Ramsey numbers and adiabatic quantum computing 
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The graph-theoretic Ramsey numbers are notoriously difficuh to calculate. In fact, for the two- 
color Ramsey numbers R{m, n) with m,n > 3, only nine are currently known. We present a quantum 
algorithm for the computation of the Ramsey numbers R{m,n). We show how the computation 
of Rljn, n) can be mapped to a combinatorial optimization problem whose solution can be found 
using adiabatic quantum evolution. We numerically simulate this adiabatic quantum algorithm and 
show that it correctly determines the Ramsey numbers ii(3, 3) and -R(2, s) for 5 < s < 7. We then 
discuss the algorithm's experimental implementation, and close by showing that Ramsey number 
computation belongs to the quantum complexity class QMA. 
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In an arbitrary party of N people one might ask 
whether there is a group of m people who are all mu- 
tually acquainted, or a group of n people who are all 
mutual strangers. Using Ramsey theory [H, Q , it can be 
shown that a threshold value R{m, n) exists for the party 
size N so that when N > R{m, n), all parties of N people 
will either contain m mutual acquaintances, or n mutual 
strangers. The threshold value R{m, n) is an example 
of a two-color Ramsey number. Other types of Ramsey 
numbers exist, though we will focus on two color Ramsey 
numbers in this paper. 

One can represent the A'^-person party problem by an 
A^-vertex graph. Here each person is associated with a 
vertex, and an edge is drawn between a pair of vertices 
only when the corresponding people know each other. In 
the case where m people are mutual acquaintances, there 
will be an edge connecting any pair of the m correspond- 
ing vertices. Similarly, if n people are mutual strangers, 
there will be no edge between any of the n correspond- 
ing vertices. In the language of graph theory Q, the m 
vertices form an m-clique, and the n vertices form an n- 
independent set. The party problem is now a statement 
in graph theory: if N > R{m, n), every graph with A'^ ver- 
tices will contain either an m-clique, or an n-independent 
set. Ramsey numbers can also be introduced using color- 
ings of complete graphs, and R{m, n) corresponds to the 
case where only two colors are used. 

Ramsey theory has found applications in mathematics, 
information theory, and theoretical computer science Q . 
An application of fundamental significance appears in the 
Paris-Harrington (PH) theorem of mathematical logic Q 
which established that a particular statement in Ramsey 
theory related to graph colorings and natural numbers 
is true, though unprovable within the axioms of Peano 
arithmetic. Such statements are known to exist as a con- 
sequence of Godel's incompleteness theorem, though the 
PH theorem provided the first natural example. Deep 
connections have also been shown to exist between Ram- 
sey theory, topological dynamics, and ergodic theory [^. 



Ramsey numbers grow extremely quickly and so are 
notoriously difficult to calculate. In fact, for two color 
Ramsey numbers R{m, n) with m,n > 3, only nine are 

presently known[3]. To check whether A'^ = R{m,ri) re- 
quires examining all 2^(^~i)/2 A^-vertex graphs. The 
number of graphs to be checked thus grows super- 
exponentially with A^, and so the task quickly becomes 
intractable. Ketonen and Solovay 0] have shown that 
this is the root cause for why the statement in the PH 
theorem cannot be proved within Peano arithmetic. 

In this paper we: (i) present a quantum algorithm for 
calculating Ramsey numbers based on adiabatic quan- 
tum evolution; (ii) numerically simulate the algorithm 
to verify that it correctly calculates small Ramsey num- 
bers; (iii) discuss its experimental implementation; and 
(iv) show that Ramsey number computation belongs to 
the quantum complexity class QMA. 

Optimization Problem: We begin by establishing a 
1-1 correspondence between the set of A^- vertex graphs 
and binary strings of length L = N{N — l)/2. To each A^- 
vertex graph G there corresponds a unique adjacency ma- 
trix A{G) which is an A'^ X A^ symmetric matrix with van- 
ishing diagonal matrix elements, and with off-diagonal 
element a^j = 1 (0) when distinct vertices i and j are 
(are not) joined by an edge. It follows that A{G) is de- 
termined by its lower triangular part. By concatenating 
column-wise the matrix elements aij appearing below 
the principal diagonal, we can construct a unique binary 
string g{G) of length L for each graph G: 

g{G) = 02,1 • • • ajv,i 0-3,2 ■ ■ ■ ajv,2 ■ • ■ ajv,jv-i- (1) 

Given the string g{G), the following procedure deter- 
mines the number of m-cliques in G. Choose m vertices 
Sa = {wi, . . . , Vm} from the A^ vertices of G and form 
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the product Ca 



) O-Vj ,Vk ■ 



Note that Ca 



when Sa forms an m-clique; otherwise Ca — 0. Now re- 
peat this procedure for all p — G{N, m) ways of choosing 
m vertices from TV vertices, and form the sum C(G) = 



2 



Sa=i^«- -^y construction, C{G) equals the number of 
m-cliques contained in G. A similar procedure deter- 
mines the number of n-independent sets in G. Briefly, 
choose n vertices Tq, = {ui,...,u„} from the N ver- 
tices in G, and form the product — Illi^^j^gT ) ''^■uj,vk j 
where ay- ^y^, = 1 — cLvj.vk ■ If the vertex set Ta forms an n- 
independent set, then = 1; otherwise = 0. Repeat 
this for all v = C{N, n) ways of choosing n vertices from 
N vertices, then form the sum I{G) — ^a- By con- 

struction, X(G) gives the number of n-independent sets 
contained in G. Finally, define 

hiG) ^C{G)+I{G). (2) 

It follows from the above discussion that h{G) is the total 
number of m-cliques and n-independent sets in G. Thus 
h{G) > for aU graphs G; and h{G) = if and only if G 
does not contain an m-clique or n-independent set. 

We can use h{G) as the cost function for the follow- 
ing combinatorial optimization problem. For given in- 
tegers (A^, m,n), and with h{G) defined as above, find 
an A'^- vertex graph G* that yields the global minimum of 
h{G). Notice that if < R{m, n), the (global) minimum 
will be /i(G,) = since Ramsey theory guarantees that 
a graph exists which has no m-clique or n-independent 
set. On the other hand, if iV > R{m,n), Ramsey theory 
guarantees h{G^,) > 0. If we begin with N < R{m,n) 
and increment iV by 1 until we first find h(G<,) > 0, then 
the corresponding N will be exactly R(m,n). We now 
show how this combinatorial optimization problem can 
be solved using adiabatic quantum evolution. 

Quantum Algorithm: The adiabatic quantum evo- 
lution (AQE) algorithm [sf exploits the adiabatic dynam- 
ics of a quantum system to solve combinatorial optimiza- 
tion problems. The AQE algorithm uses the optimization 
problem cost function to define a problem Hamiltonian 
Hp whose ground-state subspace encodes all problem so- 
lutions. The algorithm evolves the state of an L-qubit 
register from the ground-state of an initial Hamiltonian 
Hi to the ground-state of Hp with probability approach- 
ing 1 in the adiabatic limit. An appropriate measurement 
at the end of the adiabatic evolution yields a solution of 
the optimization problem almost certainly. The time- 
dependent Hamiltonian H{t) for global AQE is 

H{t) ^(^l-l^H,+ (^l^ Hp, (3) 

where T is the algorithm runtime, and adiabatic dynam- 
ics corresponds to T cx). 

To map the optimization problem associated with com- 
puting R{m, n) onto an adiabatic quantum computation, 
we begin with the 1-1 correspondence between iV- vertex 
graphs G and length L = N{N — l)/2 binary strings 
g{G). From Eq. ([T]) we see that position along the string 
is indexed by vertex pairs We thus identify a qubit 

with each such pair and will thus need L qubits. 



Defining the computational basis states (CBS) to be the 
eigenstates of cr^ (g) • • • (g) crjf~^, we identify the 2^ graph 
strings g{G) with the 2^ CBS: g{G) \g{G)). The 
problem Hamiltonian Hp is defined to be diagonal in 
the computational basis with eigenvalue h(G) associated 
with eigenstate \g{G)): 

Hp\g{G))^h{G)\g{G)). (4) 

Note that the ground-state energy of Hp will be zero iff 
there is a graph with no m-cliques or n-independent sets. 
We give an operator expression for Hp below. The initial 
Hamiltonian Hi is chosen to be 

H^ = Y.\{l'-<)^ (5) 

1=0 

where /' and cr^ are the identity and x-Pauli operator for 
qubit I, respectively. The ground-state of Hi is the easily 
constructed uniform superposition of CBS. 

The quantum algorithm for computing R(m, n) begins 
by setting TV equal to a strict lower bound for R{m, n) 
which can be found using the probabilistic method or 
a table of two-color Ramsey numbers . The AQE al- 
gorithm is run on = N{N — l)/2 qubits, and the 
energy E is measured at the end of algorithm execu- 
tion. In the adiabatic limit the result will he E — 
since N < R{m,n). The value of A'' is now incremented 
A^ — > A^+1, the AQE algorithm is re-run on L^+i qubits, 
and the energy E measured at the end of algorithm ex- 
ecution. This process is repeated until E > first oc- 
curs, at which point the associated A' will be equal to 
R{m,n). Note that any real application of AQE will 
only be approximately adiabatic. Thus the probability 
that the measured energy E will be the ground-state en- 
ergy will be 1 — e. In this case, the algorithm must be 
run k ^ 0(ln[l — 5]/lne) times so that, with probabil- 
ity (5 > 1 — e, at least one of the measurement outcomes 
will be the true ground-state energy. We can make S 
arbitrarily close to 1 by choosing k sufficiently large. 

Simulation Results: To test the adiabatic quantum 
computation of R{m,n), we numerically simulated the 
Schrodinger dynamics generated by the AQE Hamilto- 
nian H{t). Clearly, these simulations can only be run at 
finite values of T. As in Ref. [13], we chose T so that 
the algorithm success probability Pg is large compared 
to the probability that a randomly chosen CBS will be- 
long to the Z3-degenerate ground-state eigenspace of Hp 
{Ps > D/2^)). Here Ps is the probability that an en- 
ergy measurement done at the final time T will yield the 
ground-state energy Egs of Hp. Since a classical com- 
puter cannot efficiently simulate the dynamics of a quan- 
tum system, we can only obtain small Ramsey numbers. 
In this case. Hp can be found by evaluating the cost func- 
tion h{G) using the procedure described above Eq. ([2]). 

We simulated the AQE computation of i?(3, 3) and 
i?(2, s) for 5 < s < 7. Straightforward arguments ^Sj] give 
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i?(3, 3) = 6 and i?(2, s) = s. We present our simulation 
results in TableHl We see that for all m, n considered, the 
threshold value where Egg > first occurs is precisely 
at the Ramsey number: Nt = R{m, n). 

For R{2, s) and N = s, TableUgives Eg^ = 1. For these 
cases, graphs corresponding to ground-states of Hp will 
thus contain either a single s-independent set or a sin- 
gle 2-clique. There is only one s-vertex graph with an 
s-independent set, and there are C(s, 2) — s{s ~ l)/2 
graphs with one 2-clique (viz. edge). Thus the ground- 
state degeneracy D = 1 + C{s, 2), in agreement with the 
R{2,s) degeneracies in Table H] for N = s = 5,6,7. For 
i?(3, 3) and iV = 6, TableUgives Eg^ ^ 2. Thus graphs 
corresponding to ground-states are those with: (i) two 3- 
cliques; (ii) two 3-independent sets; or (iii) one 3-clique 



and one 3-independent set. Ref. derived the min- 
imum number of 3-cliques and 3-independent sets that 
can be present in an A^-vertex graph. This minimum is 
precisely our Egg for i?(3, 3) and a given A^. For = 6, 
the minimum value is 2, in agreement with Egs = 2 in Ta- 
ble m We carried out both analytical [l^ and numerical 
counts of the ground-state graphs for R{3, 3) and iV = 6. 
Both approaches found 1760 graphs giving a ground-state 
degeneracy D — 1760. In all cases appearing in Table U 
the upward jump in D seen upon reaching the Ramsey 
threshold N — R{m,n) (from below) is responsible for 
the jump in the success probability Ps also seen at this 
threshold. 

Although we would like to have calculated larger Ram- 
sey numbers, this was simply not practical. Note that 
the N ^ 7 simulations use L = 21 qubits. These sim- 
ulations are at the upper limit of 20-22 qubits at which 
simulation of the full AQE Schrodinger dynamics is prac- 
tical [13, [3 14 1 . The next smallest Ramsey number is 



i?(2,8) = 8 which requires a 28 qubit simulation, well 
beyond what can be done practically. 

Experimental Implementation: We begin by de- 
termining an operator expression for the problem Hamil- 
tonian Hp which then fixes the AQE Hamiltonian H{t) 
through Eqs. ([3]) and ([5]). Recall that the eigenvalue 
h{G) = C{G) + I{G) counts the total number of m- 
cliques and n-independent sets in a graph G. For an 
TO-vertex set Sa — {vi, ■ • ■ , Wm}, we define the edge set 
Ea = {e^ : k = 1, • . • , C(m, 2)} as the set of all edges 
connecting pairs of vertices Vi,Vj S Sa, and C(m, 2) is the 
number of ways of choosing 2 vertices out of m. If Sa cor- 
responds to an TO-clique in the graph G, the graph-string 
g{G) must have I's at all bit-positions associated with the 
edges of Ea- Let the states |0) and |1) satisfy az\a) — 
{—l)°-'\a). Then the operator Ha = nee_B (where 
Pf = (1/2) [I'^' — erf], and e labels the qubit associated 
with edge e) will have \g{G)) as an eigenstate with eigen- 
value 1 when Sa is an m-clique, and zero otherwise. The 
operator that counts all m-cliques in a graph G is then 
^cT = EaL't'^'^Ha, and by construction, H^\giG)) = 
C{G)\g{G)). A similar analysis can be carried out for n- 



independent sets. Let Ta = {vi, . . . , Vn} be an arbitrary 
n- vertex set, and Ea its corresponding edge set. If Ta 
is an n-independent set in a graph G, then the graph- 
string g{G) must have O's at all bit-positions associated 
with the edges of Ea- The operator Ha — Ileei? 
(where Pq = (1/2) [I'^ + erf], and e labels the qubit as- 
sociated with edge e) will have eigenstate \g{G)) with 
eigenvalue 1 (0) when Ta is (is not) an n-independent 
set. The operator that counts all n-independent sets in 

an arbitrary graph G is then H"^ = X^aiT'"'' ^<^' ^^"^ 
construction, H"g\g{G)) — I{G)\g{G)). For calculation 
of R{m,n), the problem Hamiltonian Hp™^ is then 



(6) 



Note that iJjy™" 



contains 0{N'') terms, where A^ is 
the number of vertices and s = max{G(A'^, m), C{N, n)}. 
Since each Ha and Ha is a projection operator, their 
operator norm will be unity and their matrix elements, 
being O's and I's, are specified with a single bit. Lastly, 
note that each term in iJp is a product of at most 
t = max{G(m, 2), G(n, 2)1 cr^-operators so that ifj^™" 
is a t-local Hamiltonian [15[ . By using perturbative gad- 



gets, it can be reduced to a 2-local Hamiltonian [16l - ll8l ]. 

For a given Hamiltonian H{t), two approaches have 
been demonstrated to experimentally implement AQE 
H-ill. Refs. [11], partitioned the fuh evolution into 
TV subintervals of duration Ai = T/M which are suffi- 
cently short that the propagator Ui for each subinterval 
I can be factored via a Trotter expansion. This approach 
was applied to three-qubit systems, though it can be used 
for arbitrary size qubit systems. Ref. |21j describes ex- 
periments using a quantum annealing device designed to 
implement adiabatic quantum optimization algorithms. 
Results are reported of AQE solution for the groundstate 
of randomly generated instances of an 8-qubit quantum 
Ising spin glass. Work using perturbative gadgets is un- 
derway to convert ifp ™" into a 2-local form amenable 
to both AQE experimental approaches. 

Ramsey Numbers and QMA: Quantum complex- 
ity theory formalizes the notion of efficient quantum algo- 
rithms. Our interest is in the quantum complexity class 
QMA which generalizes the randomized version of the 
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classical complexity class NP isl . 

QMA is a class of promise problems where each prob- 
lem L is the union of two disjoint sets of binary strings 
Ly and L„ corresponding to Yes and No instances of the 
problem. For a string x G Ly U Ln, the task is to de- 
termine whether x G Ly or x ^ Ln using polynomial re- 
sources. Let H denote a two-dimensional Hilbert space; 
and \x) the CBS labeled by the binary string x. 

Definition 1 (QMA) Let x E L = Ly U Ln and 

e = 2^^*^l^l). The promise problem L belongs to 
QMA if there exists a quantum polynomial-time veri- 
fier V{\x),\y)) {0,1}, and a polynomial 7r(|a:|) such 
that: (i) for all x € Ly, there exists an |^) G H'^d^D 
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TABLE I: Simulation results for Ramsey numbers R{3, 3) and R(2, s) for 5 < s < 7. Here is the number of graph vertices; 
Egs and D are the ground-state energy and degeneracy, respectively, for the problem Hamiltonian Hp; and T and Ps are, 
respectively, the algorithm runtime and success probability. 
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such that Pr{V{\x), |^)) = 1} > 1 — e; and (ii) for all 
X & Ln and 10 e 'H'^^l'^l^ Pr{y(|a;), |0) = !}<£. 77ere 
Pr{y(|a:), |0) = 1} is the probability that V concludes 
X G Ly when the quantum witness is \^). 

Informally, if a; is a Yes (No) instance, there exists a 
(no) quantum witness which causes V to correctly 
(mistakenly) conclude x ^ Ly with probability at least 
1 — e (greater than e) . 

A promise problem is (5AL4-Complete if it belongs to 
QMA and all problems in QMA are polynomially re- 
ducible to it. It has been shown [l^, [iSi] that k-Local 
Hamiltonian is QMA-Complete for k > 2. 

Definition 2 (fc-Local Hamiltonian) Consider an L- 
qubit Hamiltonian H = X]j=i^j' where r — poly{L); 
and each term Hj acts on at most k qubits (k-local); has 
operator norm \\Hj\\ < poly{L); and matrix elements 
specified by poly{L) bits. Finally, two constants a < h 
are specifed. The Hamiltonian H is a Yes instance if 
its groundstate energy Egs < a, and a No instance if 
Egs > b. The problem is, given a k-local Hamiltonian H , 
determine whether H is a Yes or a No instance. 

Our Ramsey number AQE algorithm leads naturally to 
an example of t-Local Hamiltonian which we call RAM- 
SEY. We have seen that the Ramsey problem Hamil- 
tonian Hp™" is a i-local Hamiltonian; is a sum of a 
polynomial number of terms Hj = Ha or Ha', and each 
Hj satisfies the polynomial bounds specified in Defini- 
tion [21 Suitable choices for the constants a and b are 
0.01 < a < 0.1 and = 1 - a. Yes instances of RAMSEY 
then correspond to < R{m, n) since Egg — Q < a, and 
No instances to A^ > R{m, n) where Egs > 1 > b. It is 
possible to carry over the proof that fc-Local Hamiltonian 
is in QMA ^ to show that RAMSEY is also in QMA. 

For an AQE algorithm with non- degenerate ground- 
state (GS), the runtime is largely determined Q by 
the minimum energy gap A — mint{£^i(i) — ii^o(i)}- 
This connection fails for the Ramsey algorithm when 
A^ — Rim, n) as the GS becomes degenerate during its 
execution and so A vanishes. Determining how the run- 
time scales when A = (as with the Ramsey algorithm) 
is an open problem in adiabatic quantum computing. 

In this paper we have presented a quantum algorithm 
that calculates two-color Ramsey numbers R{m,n); nu- 



merically simulated the algorithm and shown it correctly 
determined small Ramsey numbers; discussed its experi- 
mental implementation; and shown that Ramsey number 
computation is in the quantum complexity class QMA. 
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